Fine mapping and candidate gene analysis of CRA3.7 conferring clubroot resistance in Brassica rapa.

KEY MESSAGE: In this study, we fine-mapped a clubroot resistance gene CRA3.7 in Chinese cabbage and developed its closely linked marker syau-InDel3008 for marker-assisted selection in CR cultivars breeding. Chinese cabbage is an important leafy vegetable rich in many nutrients widely grown in China. Clubroot disease caused by an obligate biotrophic pathogen Plasmodiophora brassicae was rapidly spread and challenged to Chinese cabbage production. A clubroot resistance (CR) gene, CRA3.7, was mapped on chromosome A03 of Brassica rapa. A Chinese cabbage line 'CR510', which harbor homozygous resistance locus CRA3.7 was selected from a BC4F3 family. 'CR510' was crossed with a clubroot susceptible Chinese cabbage inbred line '59-1'. Total 51 recombinant plants were identified from an F2 population including 3000 individuals. These recombinants were selfed and the clubroot resistance of F2/3 families was evaluated. Finally, a clubroot resistance gene CRA3.7 was fine-mapped to an interval of approximately 386 kb between marker syau-InDel3024 and syau-InDel3008. According to the reference genome, total 54 genes including five encoding the TIR-NBS-LRR proteins was annotated in the fine-mapped region. Further, nine candidate's gene expression in parental lines at 7, 14 and 21 days after inoculation of P. brassicae were evaluated. Bra019376, Bra019401, Bra019403 and Bra019410 are highly expressed in 'CR510' than '59-1'. Gene sequence of Bra019410 from 'CR510' was cloned and identified different from CRa. Therefore, Bra019376, Bra019401, Bra019403 and Bra019410 are the most likely candidates for CRA3.7. Our research provides a valuable germplasm resource against P. brassicae Pb3 and CRA3.7 closely linked marker for marker-assisted selection in CR cultivars breeding.

Identification of Novel Locus RsCr6 Related to Clubroot Resistance in Radish (Raphanus sativus L.).

Clubroot is a devastating disease that causes substantial yield loss worldwide. However, the inheritance and molecular mechanisms of clubroot resistance during pathogen infection in radish remain largely unclear. In this study, we investigated the inheritance of clubroot resistance in the F2 population derived from crossing clubroot-resistant (CR) and clubroot-susceptible inbred lines "GLX" and "XNQ," respectively. Genetic analysis revealed that a single dominant gene controlled the clubroot resistance of "GLX" with a Mendelian ratio of resistance and susceptibility of nearly 3:1. Bulked segregant analysis combined with whole-genome resequencing (BSA-seq) was performed to detect the target region of RsCr6 on chromosome Rs8. Linkage analysis revealed that the RsCr6 locus was located between two markers, HB321 and HB331, with an interval of approximately 92 kb. Based on the outcomes of transcriptome analysis, in the RsCr6 locus, the R120263140 and R120263070 genes with a possible relation to clubroot resistance were considered candidate genes. In addition, three core breeding materials containing the two reported quantitative trait loci (QTLs) and our novel locus RsCr6 targeting clubroot resistance were obtained using marker-assisted selection (MAS) technology. This study reveals a novel locus responsible for clubroot resistance in radishes. Further analysis of new genes may reveal the molecular mechanisms underlying the clubroot resistance of plants and provide a theoretical basis for radish resistance breeding.

Starch content changes and metabolism-related gene regulation of Chinese cabbage synergistically induced by Plasmodiophora brassicae infection.

Clubroot is one of the major diseases adversely affecting Chinese cabbage (Brassica rapa) yield and quality. To precisely characterize the Plasmodiophora brassicae infection on Chinese cabbage, we developed a dual fluorescent staining method for simultaneously examining the pathogen, cell structures, and starch grains. The number of starch (amylopectin) grains increased in B. rapa roots infected by P. brassicae, especially from 14 to 21 days after inoculation. Therefore, the expression levels of 38 core starch metabolism genes were investigated by quantitative real-time PCR. Most genes related to starch synthesis were up-regulated at seven days after the P. brassicae inoculation, whereas the expression levels of the starch degradation-related genes increased at 14 days after the inoculation. Then genes encoding the core enzymes involved in starch metabolism were investigated by assessing their chromosomal distributions, structures, duplication events, and synteny among Brassica species. Genome comparisons indicated that 38 non-redundant genes belonging to six core gene families related to starch metabolism are highly conserved among Arabidopsis thaliana, B. rapa, Brassica nigra, and Brassica oleracea. Genome sequencing projects have revealed that P. brassicae obtained host nutrients by manipulating plant metabolism. Starch may serve as a carbon source for P. brassicae colonization as indicated by the histological observation and transcriptomic analysis. Results of this study may elucidate the evolution and expression of core starch metabolism genes and provide researchers with novel insights into the pathogenesis of clubroot in B. rapa.

A Loop-Mediated Isothermal DNA Amplification (LAMP) Assay for Detection of the Clubroot Pathogen Plasmodiophora brassicae.

Clubroot caused by Plasmodiophora brassicae is a serious threat to cruciferous crops around the world. The resting spores of P. brassicae are a primary source of infection and can survive in soil for many years. Detection of resting spores in soil is essential for forecasting clubroot prevalence. Detection of P. brassicae has been relying on plant bioassays or PCR-based methods. The loop-mediated isothermal DNA amplification (LAMP) is a promising approach for microorganism detection with the advantage of high sensitivity, accuracy, and convenience in viewing. In this study, we developed a LAMP assay for detection of P. brassicae in soil, roots, and seeds. This method can detect P. brassicae at a minimal amount of 1 fg of plasmid DNA or 10 resting spores in the soil. Compared with conventional PCR, the LAMP was more sensitive in detection of P. brassicae at the lower levels in soil samples. In conclusion, we elaborated a sensitive, accurate, and easy-to-use LAMP assay to detect P. brassicae, which will facilitate sustainable clubroot management and planning.

Spatiotemporal Quantification of Plasmodiophora brassicae Inoculum in Relation to Clubroot Development Under Inoculated and Naturally Infested Field Conditions.

Clubroot caused by Plasmodiophora brassicae is a destructive disease of cruciferous plants worldwide. A quantitative PCR (qPCR) system specific to P. brassicae was developed. Analysis of the qPCR sensitivity indicated that the lower limit of detection was 1 × 101 resting spores/ml, 1 × 102 spores/g of soil, and 1 × 103 spores/g of roots and seeds. The regression curves generated from the qPCR data of different samples had a parallel relationship. The difference between the theoretical and actual concentrations was lowest at 1 × 105 spores/g of sample, compared with other concentrations. The P. brassicae biomass in soil and plant root tissues after inoculated with different spore concentrations was correlated. A correlation analysis confirmed that the clubroot incidence and disease index at 6 weeks after inoculation increased as the spore concentration increased. Under field conditions, the natural inoculum density of the P. brassicae population decreased at the early stage and then increased, with P. brassicae mainly being detected at a soil depth of 0 to 50 cm. The horizontal distribution of P. brassicae varied in the field with occurrences of hot spots. This study established a qPCR-based method for quantitative detection of clubroot. The developed assay is useful for monitoring the spatiotemporal dynamics of P. brassicae in the field. It may also be applicable for clubroot forecasting as a part of proactive disease management.

Development of a Sinitic Clubroot Differential Set for the Pathotype Classification of Plasmodiophora brassicae.

Plasmodiophora brassicae, which is known for its broad genetic diversity for virulence, is the causal agent of clubroot disease of Brassica crops worldwide. Studies on pathotype characterization with four differential hosts according to Williams' classification system showed the predominance of pathotype 4 in China. However, the genetic variability within pathotype 4 complicates the breeding of durable clubroot-resistant (CR) cultivars. Herein, a Sinitic clubroot differential (SCD) set was developed using a set of eight differential inbred lines of Chinese cabbage with known or novel CR genes. The presence of immense diversity within pathotype 4 of Williams' system was verified, and 11 pathotypes were characterized using the developed SCD system. The scalability and practicability of the system was further confirmed with a subset of 95 field isolates from different Brassica crops and different regions of China and Korea. Sixteen pathotypes were detected from 132 field isolates, named Pb1 to Pb16, respectively. Among them, Pb1 and Pb4 were prevalent in diverse Brassica crops in the southern and northern regions of China. Pb12, Pb13, Pb14, and Pb16 showed area-specific distribution. The SCD set developed herein will provide important genetic resources for pathogenicity studies of P. brassicae and for CR breeding in Chinese cabbage and other Brassica crops.

Mining of Brassica-Specific Genes (BSGs) and Their Induction in Different Developmental Stages and under Plasmodiophora brassicae Stress in Brassica rapa.

Orphan genes, also called lineage-specific genes (LSGs), are important for responses to biotic and abiotic stresses, and are associated with lineage-specific structures and biological functions. To date, there have been no studies investigating gene number, gene features, or gene expression patterns of orphan genes in Brassica rapa. In this study, 1540 Brassica-specific genes (BSGs) and 1824 Cruciferae-specific genes (CSGs) were identified based on the genome of Brassica rapa. The genic features analysis indicated that BSGs and CSGs possessed a lower percentage of multi-exon genes, higher GC content, and shorter gene length than evolutionary-conserved genes (ECGs). In addition, five types of BSGs were obtained and 145 out of 529 real A subgenome-specific BSGs were verified by PCR in 51 species. In silico and semi-qPCR, gene expression analysis of BSGs suggested that BSGs are expressed in various tissue and can be induced by Plasmodiophora brassicae. Moreover, an A/C subgenome-specific BSG, BSGs1, was specifically expressed during the heading stage, indicating that the gene might be associated with leafy head formation. Our results provide valuable biological information for studying the molecular function of BSGs for Brassica-specific phenotypes and biotic stress in B. rapa.

Identification and Mapping of the Clubroot Resistance Gene CRd in Chinese Cabbage (Brassica rapa ssp. pekinensis).

The rapid spread of clubroot disease, which is caused by Plasmodiophora brassicae, threatens Brassicaceae crop production worldwide. Breeding plants that have broad-spectrum disease resistance is one of the best ways to prevent clubroot. In the present study, eight Chinese cabbage germplasms were screened using published clubroot-resistant (CR) loci-/gene-linked markers. A CR gene Crr3 potential carrier "85-74" was detected which linked to marker BRSTS61; however, "85-74" shows different responses to local pathogens "LAB-19," "LNND-2," and "LAB-10" from "CR-73" which harbors Crr3. We used a next-generation sequencing-based bulked segregant analysis approach combined with genetic mapping to detect CR genes in an F2 segregant population generated from a cross between the Chinese cabbage inbred lines "85-74" (CR) and "BJN3-1" (clubroot susceptible). The "85-74" line showed resistance to a local pathogen "LAB-19" which was identified as race 4; a genetic analysis revealed that the resistance was conferred by a single dominant gene. The CR gene which we named CRd was mapped to a 60 kb (1 cM) region between markers yau389 and yau376 on chromosome A03. CRd is located upstream of Crr3 which was confirmed based on the physical positions of Crr3 linked markers. The identification of CRd linked markers can be applied to marker-assisted selection in the breeding of new CR cultivars of Chinese cabbage and other Brassica crops.

Integrated analysis of leaf morphological and color traits in different populations of Chinese cabbage (Brassica rapa ssp. pekinensis).

KEY MESSAGE: QTLs and candidate gene markers associated with leaf morphological and color traits were identified in two immortalized populations of Brassica rapa, which will provide genetic information for marker-assisted breeding. Brassica rapa is an important leafy vegetable consumed worldwide and morphology is a key character for its breeding. To enhance genetic control, quantitative trait loci (QTLs) for leaf color and plant architecture were identified using two immortalized populations with replications of 2 and 4 years. Overall, 158 and 80 QTLs associated with 23 and 14 traits were detected in the DH and RIL populations, respectively. Among them, 23 common robust-QTLs belonging to 12 traits were detected in common loci over the replications. Through comparative analysis, five crucifer genetic blocks corresponding to morphology trait (R, J&U, F and E) and color trait (F, E) were identified in three major linkage groups (A2, A3 and A7). These might be key conserved genomic regions involved with the respective traits. Through synteny analysis with Arabidopsis, 64 candidate genes involved in chlorophyll biosynthesis, cell proliferation and elongation were co-localized within QTL intervals. Among them, SCO3, ABI3, FLU, HCF153, HEMB1, CAB3 were mapped within QTLs for leaf color; and CYCD3;1, CYCB2;4, AN3, ULT1 and ANT were co-localized in QTL regions for leaf size. These robust QTLs and their candidate genes provide useful information for further research into leaf architecture with crop breeding.

Genome-Wide Analysis and Characterization of Aux/IAA Family Genes in Brassica rapa.

Auxins are the key players in plant growth development involving leaf formation, phototropism, root, fruit and embryo development. Auxin/Indole-3-Acetic Acid (Aux/IAA) are early auxin response genes noted as transcriptional repressors in plant auxin signaling. However, many studies focus on Aux/ARF gene families and much less is known about the Aux/IAA gene family in Brassica rapa (B. rapa). Here we performed a comprehensive genome-wide analysis and identified 55 Aux/IAA genes in B. rapa using four conserved motifs of Aux/IAA family (PF02309). Chromosomal mapping of the B. rapa Aux/IAA (BrIAA) genes facilitated understanding cluster rearrangement of the crucifer building blocks in the genome. Phylogenetic analysis of BrIAA with Arabidopsis thaliana, Oryza sativa and Zea mays identified 51 sister pairs including 15 same species (BrIAA-BrIAA) and 36 cross species (BrIAA-AtIAA) IAA genes. Among the 55 BrIAA genes, expression of 43 and 45 genes were verified using Genebank B. rapa ESTs and in home developed microarray data from mature leaves of Chiifu and RcBr lines. Despite their huge morphological difference, tissue specific expression analysis of BrIAA genes between the parental lines Chiifu and RcBr showed that the genes followed a similar pattern of expression during leaf development and a different pattern during bud, flower and siliqua development stages. The response of the BrIAA genes to abiotic and auxin stress at different time intervals revealed their involvement in stress response. Single Nucleotide Polymorphisms between IAA genes of reference genome Chiifu and RcBr were focused and identified. Our study examines the scope of conservation and divergence of Aux/IAA genes and their structures in B. rapa. Analyzing the expression and structural variation between two parental lines will significantly contribute to functional genomics of Brassica crops and we belive our study would provide a foundation in understanding the Aux/IAA genes in B. rapa.

Quantitative Trait Loci for Morphological Traits and their Association with Functional Genes in Raphanus sativus.

Identification of quantitative trait loci (QTLs) governing morphologically important traits enables to comprehend their potential genetic mechanisms in the genetic breeding program. In this study, we used 210 F2 populations derived from a cross between two radish inbred lines (Raphanus sativus) "835" and "B2," including 258 SSR markers were used to detect QTLs for 11 morphological traits that related to whole plant, leaf, and root yield in 3 years of replicated field test. Total 55 QTLs were detected which were distributed on each linkage group of the Raphanus genome. Individual QTLs accounted for 2.69-12.6 of the LOD value, and 0.82-16.25% of phenotypic variation. Several genomic regions have multiple traits that clustered together, suggested the existence of pleiotropy linkage. Synteny analysis of the QTL regions with A. thaliana genome selected orthologous genes in radish. InDels and SNPs in the parental lines were detected in those regions by Illumina genome sequence. Five identified candidate gene-based markers were validated by co-mapping with underlying QTLs affecting different traits. Semi-quantitative reverse transcriptase PCR analysis showed the different expression levels of these five genes in parental lines. In addition, comparative QTL analysis with B. rapa revealed six common QTL regions and four key major evolutionarily conserved crucifer blocks (J, U, R, and W) harboring QTL for morphological traits. The QTL positions identified in this study will provide a valuable resource for identifying more functional genes when whole radish genome sequence is released. Candidate genes identified in this study that co-localized in QTL regions are expected to facilitate in radish breeding programs.

Anatomic Characteristics Associated with Head Splitting in Cabbage (Brassica oleracea var. capitata L.).

Cabbage belonging to Brassicaceae family is one of the most important vegetables cultivated worldwide. The economically important part of cabbage crop is head, formed by leaves which may be of splitting and non-splitting types. Cabbage varieties showing head splitting causes huge loss to the farmers and therefore finding the molecular and structural basis of splitting types would be helpful to breeders. To determine which anatomical characteristics were related to head-splitting in cabbage, we analyzed two contrasting cabbage lines and their offspring using a field emission scanning electron microscope. The inbred line "747" is an early head-splitting type, while the inbred line "748" is a head-splitting-resistant type. The petiole cells of "747" seems to be larger than those of "748" at maturity; however, there was no significant difference in petiole cell size at both pre-heading and maturity stages. The lower epidermis cells of "747" were larger than those of "748" at the pre-heading and maturity stages. "747" had thinner epidermis cell wall than "748" at maturity stage, however, there was no difference of the epidermis cell wall thickness in the two lines at the pre-heading stage. The head-splitting plants in the F1 and F2 population inherited the larger cell size and thinner cell walls of epidermis cells in the petiole. In the petiole cell walls of "747" and the F1 and F2 plants that formed splitting heads, the cellulose microfibrils were loose and had separated from each other. These findings verified that anomalous cellulose microfibrils, larger cell size and thinner-walled epidermis cells are important genetic factors that make cabbage heads prone to splitting.

Quantitative trait loci mapping of partial resistance to Diamondback moth in cabbage (Brassica oleracea L).

KEY MESSAGE: The resistance to Diamondback moth insect in cabbage is governed by many minor loci in quantitative nature, and at least four genetic loci should be incorporated in marker-assisted breeding program for developing partially resistant DBM cabbage cultivars. The Diamondback moth (DBM), Plutella xylostella (L.), is the most destructive insect infesting cruciferous plants worldwide. Earlier studies have reported that the glossy leaves of cabbage are associated with resistance to this insect. However, until now, genetics of DBM resistance has not been studied in detail, and no QTL/gene mapping for this trait has been reported. In this paper, we report quantitative trait loci (QTL) mapping of DBM-resistant trait using 188 randomly selected segregating F 3 population derived from crossing a partially DBM-resistant glossy leaf cabbage (748) with a susceptible smooth cabbage line (747). Quantitative trait loci mapping using phenotypic data of four consecutive years (2008, 2009, 2010, and 2011) on DBM insect infestation detected a total of eight QTL on five linkage groups suggesting that DBM resistance is a quantitative in nature. Of these QTL, four QTL, i.e., qDbm 1 on LG1, qDbm5 and qDbm6 on LG7, and qDbm8 on LG9, were detected in different tests and years. The QTL, qDbm6 on LG7, was consecutively detected over 3 years. Tightly linked molecular markers have been developed for qDbm8 QTL on LG9 which could be used in marker-assisted breeding program. Our research demonstrated that for desired DBM resistance cultivar breeding, those four genetic loci have to be taken into consideration. Furthermore, the comparative study revealed that DBM resistance QTL is conserved between close relative model plant Arabidopsis thaliana and Brassica oleracea genome.

The effects of fluctuating culture temperature on stress tolerance and antioxidase expression in Esteya vermicola.

The endoparasitic nematophagous fungus, Esteya vermicola, has shown great potential as a biological control agent against the pine wood nematode, Bursaphelenchus xylophilus. Fluctuating culture temperatures can affect fungal yields and fungal tolerance to desiccation, UV radiation, H2O2, and heat stress, as well as antioxidase expression. To explore these effects, E. vermicola cultured under five temperature ranges, 26°C, 15-26°C, 26-35°C, 20-30°C, and 15-35°C, were compared. The cultures grown at lower temperatures showed better growth, stronger tolerance to desiccation, UV, and H2O2 stresses, and increased catalase expression, However, these cultures also showed weaker heat stress tolerance and lower superoxide dismutase expression than the higher-temperature cultures. In particular, the E. vermicola cultured at 20-30°C, i.e., fluctuating in a narrow range around the optimal temperature, showed the best performance. Therefore, for production in practical applications, this narrowly fluctuating, moderate temperature appears to be optimal for yield and stress tolerance in E. vermicola.

Transcriptome Analysis of Brassica rapa Near-Isogenic Lines Carrying Clubroot-Resistant and -Susceptible Alleles in Response to Plasmodiophora brassicae during Early Infection.

Although Plasmodiophora brassicae is one of the most common pathogens worldwide, the causal agent of clubroot disease in Brassica crops, resistance mechanisms to it are still only poorly understood. To study the early defense response induced by P. brassicae infection, a global transcriptome profiling of the roots of two near-isogenic lines (NILs) of clubroot-resistant (CR BJN3-2) and clubroot-susceptible (BJN3-2) Chinese cabbage (Brassica rapa) was performed by RNA-seq. Among the 42,730 unique genes mapped to the reference genome of B. rapa, 1875, and 2103 genes were found to be up- and down-regulated between CR BJN3-2 and BJN3-2, respectively, at 0, 12, 72, and 96 h after inoculation (hai). Functional annotation showed that most of the differently expressed genes are involved in metabolism, transport, signal transduction, and defense. Of the genes assigned to plant-pathogen interactions, 151 showed different expression patterns between two NILs, including genes associated with pathogen-associated molecular patterns (PAMPs) and effectors recognition, calcium ion influx, hormone signaling, pathogenesis-related (PR) genes, transcription factors, and cell wall modification. In particular, the expression level of effector receptors (resistance proteins), PR genes involved in salicylic acid (SA) signaling pathway, were higher in clubroot-resistant NIL, while half of the PAMP receptors were suppressed in CR BJN3-2. This suggests that there was a more robust effector-triggered immunity (ETI) response in CR BJN3-2 and that SA signaling was important to clubroot resistance. The dataset generated by our transcriptome profiling may prove invaluable for further exploration of the different responses to P. brassicae between clubroot-resistant and clubroot-susceptible genotypes, and it will strongly contribute to a better understanding of the molecular mechanisms of resistance genes of B. rapa against P. brassicae infection.